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Abstract — Effective sensor scheduling requires the considera- 
tion of long-term effects and thus optimization over long time 
horizons. Determining the optimal sensor schedule, however, 
is equivalent to solving a binary integer program, which is 
computationally demanding for long time horizons and many 
sensors. For linear Gaussian systems, two efficient multi-step 
sensor scheduling approaches are proposed in this paper. The 
first approach determines approximate but close to optimal 
sensor schedules via convex optimization. The second approach 
combines convex optimization with a branch-and-bound search 
for efficiently determining the optimal sensor schedule. 

I. Introduction 

The recent advances in miniaturization, wireless commu- 
nication, and sensor technology make it possible to build 
up and deploy sensor systems for a smart and persistent 
surveillance. For instance, sensor networks consisting of nu- 
merous inexpensive sensors are a popular subject in research 
and practice for monitoring physical phenomena including, 
temperature and humidity distributions, biochemical concen- 
trations, or vibrations in buildings 11]. For many of such sensor 
systems it is necessary to balance between maximizing the 
information gain and minimizing the consumption of limited 
resources like energy, computing power, or communication 
bandwidth. Sensor scheduling, which is also referred to as 
sensor selection, allows trading off these conflicting goals and 
forms the basis for an efficient and intelligent processing of 
the sensor data. 

In this paper, the sensor scheduling problem for linear 
Gaussian systems is studied, where one out of a set of 
sensors is selected at each time instant for performing a 
measurement. The main objective is to allocate the sensors 
in a most informative way, which requires making decisions 
involving multiple time steps ahead. Many of the existing 
multi-step sensor scheduling approaches for Unear Gaussian 
problems are focused on efficiently traversing a decision tree 
of sensor schedules. In order to avoid enumerating all possible 
sensor schedules of the tree, optimal or suboptimal search 
techniques are employed. While optimal techniques yield the 
optimal sensor schedule by all means (see e.g. |l2|, 13), 
suboptimal methods as those in ID allow more significant 
savings in computational demands by abdicating the guarantee 
of conserving the optimal schedule. 

Alternatively to traversing the decision tree, which cor- 
responds to solving a binary integer program, convex opti- 



mization approaches have recently been proposed for solving 
sensor selection problems, i.e., problems of selecting the 
best n-element subset from a set of sensors (see JS], IS). 
These approaches can significantly improve the efficiency of 
determining informative sensor schedules, but they are so far 
not appropriate for optimal multi-step sensor scheduling for 
arbitrary linear Gaussian dynamics and sensor models. 

Both multi-step sensor scheduling approaches proposed in 
this paper overcome these restrictions. For linear Gaussian 
systems, the sensor scheduling problem is stated in Sec. Inland 



is formulated as a binary integer problem in Sec. Ill In Sec. IV 



it is formally proven that this optimization problem is a convex 
optimization problem when employing continuous relaxation 
of the decision variables. The first approach directly solves the 
resulting convex program, which leads to suboptimal but valu- 
able sensor schedules without demanding many computations 
and memory. In order to provide the optimal sensor sequence, 
the second approach described in Sec. M utilizes branch- 
and-bound search for traversing a decision tree. To exclude 
complete subtrees containing suboptimal sensor schedules as 
early as possible, the solution of the convex optimization 
is used for calculating tight lower and upper bounds to the 
subtrees' values. The performance of the proposed approaches 



is demonstrated by means of simulations in Sec. VI while in 
Sec. I VII| conclusions and an outlook to future work are given. 

II. Problem Formulation 

In this paper, the sensor scheduling problem for discrete- 
time linear Gaussian models is examined. The dynamics model 
of the observed system is given by 
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where k ^ 0,1, . . . , is the discrete time index. A finite set S 
of sensors is considered for performing measurements, where 
measurement ^^ from sensor i e 5 = {1, . . . , S*} is related to 
the system state a;^, via the measurement model 
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Both Afc and H). are time-variant matrices. The noise terms 
w^. and vl, are zero-mean white Gaussian with covariance 
matrices C^ and C^'*, respectively. A measurement value z]. 
of sensor z € 5 is a realization of z^. The initial system state 



ajQ ^ A/'(a;Q; Xq, eg) at time step A: = is Gaussian with 
mean Xq and covariance Cg. 

The aim of muhi-step sensor scheduUng is to minimize 
the state covariance Cf, and thus the uncertainty of the 
state estimate under the consideration of the future behavior 
of the observed dynamical system and long-term sensing 
costs. For this purpose, the optimal sensor schedule u^.jy = 

[{u*if,---Au*NfV ^ {0, l}^'^ has to be determined 
over a finite iV-step time horizon. Here, the binary vector 
Mfe — [wfc,i, . . . , Wfe,s]^ encodes the index of the sensor to 
be scheduled for measurement at time step k, i.e., if sensor i 
is scheduled at time step k then Uk,i = 1 and Ukj = for all 

III. Constrained Optimization 

To determine the optimal sensor schedule m^.^ over a time 
horizon of length iV, the sensor scheduling problem for the 
problem setting given in Sec. Ill] can generally be formulated 
as constrained optimization problem according to 
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where c^, — [ck.i, ■ ■ ■ ,Ck.s]'^ contains the sensor costs c^^i, 
e.g., energy or communication, of selecting sensor i at time 
step k and J(mj.^) = J2k=i 9k{ui:k) i^ the cumulative 
objective function to be minimized. The scalar functions 5fc( • ), 
i.e., the summands of J(ui.jv)^ can be 

• the trace operator trace (C^(m]^.j,)) , 

• the root-determinant \/\Ck(—i-k)\ ' 

• or the maximum eigenvalue Amax (C^(u]^.j,)) 

of the state covariance and thus quantify the uncertainty 
subsumed in C^(uj^.j,). Due to the linearity of ([TJ and Q, the 
state covariance itself is given by the Kalman filter covariance 
recursion (given here in information form) 



j2uk,-{ui)^-{cn 
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commencing from Cg. The independence of the covariance 
recursion from the measurement values z^, i = 1,...,S, 
allows the off-line calculation of the optimal sensor schedule. 
With the contraint in Q it is guaranteed that a feasible 
sensor schedule does not exceed a maximum sensor cost C. 
The constraints in (|5]l and (|6| together ensure that one sensor 
per time step is selected for measurement. The restriction to 
one sensor per time step is made for brevity and clarity reasons 
only. The extension to selecting multiple sensors per time step 
can be achieved by replacing the right hand side of (|5]l with 
the desired number of sensors. Alternatively, by modifying 



(|3]l and Q, is also possible to minimize the sensor costs 
regarding a maximum allowed value of J( • ), i.e., a maximum 
allowed uncertainty. The results presented in the following can 
be easily altered to cover this modified optimization problem. 

IV. Convex Relaxation 

The optimization problem in ([3|-(|6]) is a so-called binary 
integer program. Problems of this type are known to be NP- 
hard (see Q) and thus, obtaining the optimal solution for large 
N and/or large S is computationally prohibitive in general. 
However, by replacing the binary non-convex constraints in 
(|6]l with the linear constraints u^. e [0, 1]'^ for fc = 1, . . . , N, 
a convex relaxation of the original problem is obtained. To 
see this, it is important to note that the constraints Q and (|5]l 
are already convex. Furthermore, as shown in the following 
theorem, the sum to be minimized in Q is now convex as 
well. 

Theorem 1 The objective function J{ui.j^) in (|3]l is convex 

in terms of Ui.^ G [0, 1]'^'^. 

Proof: To prove the convexity of gk{ui-k) ^i^^ thus of 
J{ui.pf), it must be shown that (see for example |[8|) 

5fc(A-Ui,fe + (l-A)-Ui^fc) < 

A-5fc(Mi:fe) + (l-A)-5fc(Mi,fe) (8) 

for fc = l,...,N,yuj^.j^,u^.,^ e [0, !]''■• ^, and VAe [0,1]. 

At first, it is proven by induction that the covariance 
recursion (|7| is a convex function of Uj^.^. The induction 
starts with Cf(ui). Defining M*^ := (H^) • (C^;^)^^ -H*^ 

and Pi(ui) := (Aq • Cg • A J + CJJ')"' + E. ^i,. ' M^ and 
utilizing the results in 191 on matrix convex functions, it 
follows from the matrix convexity property of the matrix 
inversion that 

Cf{X-u, + (1-A) -Ml) = (A-Pi(mi) + (1-A) -Pife))"' 
<A.Pr^fe)+(l-A).Pr^fe) 
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'^u^,u-^ e [0,1]'^ and VA G [0,1]. Defining the predicted 
covariance C^(Mi:fc_i) :== Ak-i ■ Cl_-^^{u^.k_^) ■ Al_^ + 
C™ j^, it generally holds that 

C^(A-Mi:fc + (1-A)-Mi:,) 

= (c^(A.Mi:fc_i + (l-A)-Ui^,_0"' 

s _ , 

+ ^(A-Ufc,, + (l-A)-%,)-MI, 
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for fc = 2, ...,7V, VMi:fc,Mi:fc e [O,!]'^''^, and VA G 
[0,1]. Here, (a) results from the induction hypothesis that 
^k-ii'iii-k-i) is convex in Ui.^._i, from the convexity of 
the matrix inversion, and from rearranging terms; (b) is the 
result of a repeated application of the convexity of the matrix 
inversion. 

As the trace is a linear matrix function and the root- 
determinant as well as the maximum eigenvalue are convex 
matrix functions (see for example |8|), the inequality in (|8]l 
holds if these three functions are applied on (|9]l. Thus, gk{ui-k) 
is convex and the nonnegative sum J^(Mijv) — Tlk^i 9k{y±i-k) 
is convex as well, which concludes the proof. ■ 

Theorem [T] forms one of the main contributions of this 
paper. It is important to note that the sensor scheduling 
problem formulated by (|3]l-(|7]i and its convex relaxation 
proven in Theorem [T] extends existing convex approaches 10, 
f6l in many ways. Instead of one-step time horizons, i.e., 
myoptic/greedy scheduling, arbitrarily long time horizons are 
possible. Furthermore, the dynamics model in ([TJ need not to 
be restricted to regular system matrices A^ and to system noise 
covariances C^ = 0. Especially the latter is of paramount 
importance for realistic sensor scheduling problems. Finally, 
there is no restriction to a specific scalar function gk{-) as in 
Is). Instead, various functions for evaluating the quality of a 
sensor schedule are considered here. 

A. Solving the Relaxed Problem 

The computational complexity of optimally solving the orig- 
inal binary integer program is in 0{S^). A variety of methods 
is available for efficiently solving the convex relaxation of the 
sensor scheduling problem, e.g., interior-point methods |8|. 
These methods typically require only a few tens of iterations 
for calculating the optimal solution even for large problem 
sizes, e.g., length of time horizon and number of sensors 
beyond 10. The computational complexity of one iteration is 
polynomial in the number of variables in Ui.j^, which is 5* • A^. 

However, the solution u'j^.^ of the convex problem only ap- 
proximates the optimal solution m|.„ of the original scheduling 
problem. More specifically, u\.^ is no longer binary and the 
objective function value j' :— J(u\.jv) is a lower bound of 
the optimal value J(uJ.^). The latter finding follows directly 
from the convexity of the relaxed problem and from the fact 
that the relaxed solution set [0, 1]'^'^ contains the binary set 
of the original problem. 

B. Conversion into Binary Solution 

In order to allow selecting sensors for measurement, Mi^v 
has to be converted into a binary vector by employing 
an appropriate conversion or rounding method. The value 
J" := J{ui-n) of '^he resulting (binary) sensor schedule 
m".^ has to be as close as possible to the optimal one in 
order to provide informative sensor measurements. In the 
following, two appropriate conversion methods are introduced. 
Independent of the chosen conversion method, the value J" of 
the converted sensor schedule provides an upper bound to the 
optimal value J(u5^.^). 



1) Sampling: Each component u), of u'l-jv ^^^ ^^ inter- 
preted as a discrete probability distribution over the set of 
sensor indices S. This is due to the constraint in (|5]), whereby 
the elements u\^,i = 1 , . . . , S* of u'j, are within the interval 
[0, 1] and sum up to one. Hence, a sensor i corresponding to 
an element u\ ■ with a large value can be considered as being 
more likely in the optimal sensor schedule than sensors with 
small values. 

To convert m\.^ into a feasible binary vector, for each k = 
I, . . . ,N a (single) sensor is randomly selected according to 
the distribution u'^,. For being feasible, the resulting converted 
schedule m".j, has to satisfy the cost constraint Q. Otherwise, 
the schedule is discarded. This procedure is repeated multiple 
times, where only the currently best feasible schedule, i.e., the 
schedule that satisfies (|4]) and provides the currently smallest 
objective function value J" is stored. The sampling-based 
conversion method can be terminated for example after a 
predefined number of trials or when the currently best value 
J" remains unchanged for a predefined number of trials. 

2) Swapping: A converted schedule Mijv ^'^^ ^^ improved 
by adapting the swapping method proposed in [6|. A modified 
sensor schedule is derived from u".jy by swapping a scheduled 
sensor with one of the unselected sensors for each time 
step fc. The choice of an unselected sensor at time step 
k is deterministically guided according to the probabiUties 
represented by u'^,, i.e., the sensors are selected for swapping in 
descending order of the values in w'j.. If the modified schedule 
is feasible and improves the objective function value J", it is 
used as starting point for the next swapping trial. 

In order to start the swapping method with a feasible 
schedule, the sensor schedule that selects at each time step k 
the sensor i — argmirij Ck.j with the smallest cost is chosen 
initially. The method must terminate because there is only 
a finite but very large number of swapping possibilities. To 
bound the computational demand, the number of swapping 
trials is limited by means of a predefined value. 

V. Optimal Scheduling 

Determining the optimal sensor schedule and thus, directly 
solving the binary integer program given by (|3]l-(|6]l can be 
considered as searching a decision tree with depth N and 
branching factor S. The problem here is that the optimal 
solution often can be found at an early stage, while the proof 
of its optimality requires evaluating most of the suboptimal 
sensor schedules, which is infeasible for large problem sizes. 
In this section, the previously introduced convex optimiza- 
tion approach is combined with efficient search methods 
for decision trees for early eliminating (pruning) suboptimal 
schedules. 

A. Branch-and-Bound 

A search technique common for classical decision problems 
like traveling-salesman or knapsack is branch-and-bound (BB) 
search. The basic idea of BB is to assign lower and upper 
bounds of the achievable objective function value to any 
visited node. Based on these bounds, nodes and thus complete 



Algorithm 1 Branch-and-Bound search algorithm utilizing 
convex optimization for calculating lower and upper bounds. 
The algorithm is initialized with Jnun — oo. 

1: For a given sensor schedule Uj^.^, do: 

2: if leaf node, i.e., k ^ N then 

3: '^min '^ J\Vl1:N) 

4: else 

i: U -^ % II List of sensors to expand 

6: for all sensors z e {1, . . . , S*} do 

7: // Ui.j. and Uk+i.i — 1 fixed 

8: if costi < C and Ji < J,nin then 

9: J} ^r- Solve convex optimization problem 

10: J" <— Calculate upper bound via conversion 

11: U^UU{i} 

12: end if 

13: end for 

14: U -i— sort(Z^) // Sort sensors based on j] 

15: for all sensors i £U do 

16: if Jl < Jinin and V j e W : J] < J] then 

17: Expand i II Set Ufc+i.i = 1, call Algorithm 1 

18: end if 

19: end for 

20: end if 



subtrees can be pruned under the guarantee that the pruned 
node is not part of the optimal sensor schedule. 

For a particular node that was reached during the search by 
employing the sensor schedule Uj^.^, e {0, l}'^ ' ^ , the objective 
function can be written according to 



J{U1:n) = JiUl:k) + JiUk+l:N) 



ik+l:NJ 
known unknown 



(10) 



where only the value of the first summand is evaluated and 
thus known. While the value of the second summand is not 
calculated yet, a lower and upper bound can be easily assigned 



to it by exploiting the results of Sec. IV-A and Sec. IV-B The 



value of the optimal solution lif.-^^ ■ i^r °f ^^^ convex relaxation 



^k + l:N 



for minimizing J{uk+i-N) serves as lower bound and the 



the value J, 

,1 



Jiy^l:, 



and the bounds j' 



J. 



k+lJ 

J(m'j._i_2.7v)' Jf ■— Ji + •^(^It+2iv) ^1'^ calculated, where 
Uk+i,i = 1 and Uk+1,3 = for all j ^ i. A node i is expanded 
only if following four requirements are fulfilled: (1) the cost 
constraint can be met, i.e., a feasible solution exists (line 8), 
(2) the value Ji of the node is below the value Jirdn of the 
currently best sensor schedule (Une 8), (3) the lower bound j] 
is below Jniin (line 16), and (4) the lower bound is below the 
upper bounds of all neighboring nodes j ^ i (line 16). 

Obviously, the third requirement implies the second one. 
But in order to avoid an unnecessary calculation of the lower 
and upper bound, the second requirement is checked separately 
together with the first requirement (line 8-12). To further 
accelerate the search, the remaining sensors in U are sorted 
in descending order according of their lower bounds (line 12). 
In doing so, the search is continued with the most promising 
sensor first in order to force a stronger reduction of the 
currently best value J„iin- This value is automatically reduced 
once a leaf node is reached (line 2-3). 

VI. Simulation Results 

The effectiveness of the proposed sensor scheduling 
methods is demonstrated in the following by means of a 
numerical simulation from the field of target tracking^ The 
state x^. = [ccfe, a;^, yk^Vk]'^ of the observed target comprises 
the two-dimensional position [xkjyk]"^ and the velocities 
[xkiVk]"^ in X and y direction. The system matrix and noise 
covariance matrix of Wk of the dynamics model ([T]) are 



and C^ 



2 

T 



(11) 



respectively, where I„ indicates an n x n identity matrix 
and (E) is the Kronecker matrix product. In ( [TT| ), T = 1 s 
is the sampling interval and q — 0.2 is the scalar diffusion 
strength. Mean and covariance of the initial state Xq are 
Xq ~ [0, 1, 0, 1]"^ and Cg = 10 • I4, respectively. 

A sensor network observes the target. It consists of six 
sensors with measurement matrices 



1 

k+l:N 



into a binary-valued vector Ufc+Liv H^ = H| ^ [l O] , H^ = H| = [O 1 O] , 



conversion of u 

provides an upper bound. Hence, the inequality 

JiUl:k) + -^("i+Lw) < JiUl-.w) < JiUl:k) + -'"(Mfc+1:A') 

holds for the objective function value in ( [TOJ I. 

B. Search Algorithm 

The combination of BB search with convex optimization is 
illustrated in Algorithm [T] which basically employs a depth- 
first search. For a given sensor schedule Uj^.^. it is checked 
which child nodes should be expanded, i.e., it is checked 
whether an element Uk+i,i, i ^ S of Uf^_^_l could be set to 
one or not. Therefore, for each child node i e 5 the minimum 
possible cost 



H|. = [0 1] , H^ = [0 1 0] 



noise variances Cl 

a 



v.b 



0.2, Cf^ 



C\ 



i>,3 



Cl 



0.1, 



cost, 



k 



Ck+1,'. 



N 

E 

n=k+2 



minc„ 



Cl- = 0.05, and costs c^ = [1, 2, 3, 2, 3, 2]^ for 
each k. Furthermore, it is also possible to omit a measurement. 
This option can be considered as having a seventh sensor 
with infinite noise variance. Performing no measurement is 
free of cost, i.e., Ckj = 0. Altogether, the set S comprises 
S ~ 7 sensors. The scalar functions gk{-) are set to the root- 
determinant for each k. 

For comparison, six different scheduling methods are con- 
sidered: 

'Further aspects associated to target tracking such as target detection, 
misses and false alarms, or track maintenance are omitted for simplicity. 
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Fig. 1. (a) Number of nodes in the decision tree when applying the branch-and-bound-based scheduling methods BBC (black lines), BBL (green), and 
BBZ (red) for different time horizons lengths N and for the cost functions C\ (dashed) and C2 (solid) in log-scale, (b) Objective function values J of the 
scheduling methods BBC (black, solid), CONVEX (green, dotted), GREEDY (red, dashed), and GREEDY* (blue, dash-dotted) for maximum cost function C\. 



CONVEX The approach described in Sec. IV which directly 

solves the convex optimization problem and employs the 

swapping method for conversion. 
BBC The BB approach described in Sec. |V] For determining 

the upper bounds via conversion, the swapping method 

is employed. 
BBL Like BBC but without utilizing upper bounds. 
BBZ BB search that employs no upper bounds and bounds 

the second summand in ( [TO] l from below with zero (see 

for example |10|). 
GREEDY In order to minimize the objective function J( • ), 

at each time step fc the sensor that leads to minimum 

function value g^ ( • ) and allows meeting the maximum 

cost constraint is scheduled. 
GREEDY* Greedy scheduling, where the scalar functions 

g^{-) are modified to 5fc(Mi:fc) = \/|C^(Mi:fc)| • (1 + 

cj ■ Uk) (see for example fH}, 113). 
For CONVEX and BBC, the number of swapping trials is set 
to 5 • TV. 

A. Comparison of Branch-and-Bound Methods 

In Fig. [T] (a), the search performance of the three BB 
methods is compared. For this purpose, two different maxi- 
mum costs Ci(iV) = round(1.5 • TV) and C2{N) = 3 • A^ are 
considered, which depend on the change of the time horizon 
length iV = 1, . . . , 10. The maximum cost function C2{N) al- 
lows sensor scheduling without omitting a measurement. With 
the proposed optimal scheduling method BBC, the number of 
nodes in the decision tree can be kept on a low level. Here, 
the search performance clearly benefits from the tight lower 
and upper bounds provided by the convex optimization and the 
conversion, respectively. This can be seen in particular for C2, 
where BBC only visits at most 92 nodes, while the complete 
decision tree contains X!fc=i '^^ "^ ^-^ ' ^^^ nodes. The higher 
number of visited nodes for cost function Ci compared to C2 
results from the effect that the more restrictive cost constraint 
provided by Ci leads to looser bounds. 



Without considering upper bounds for pruning as it is 
the case for BBL, the number of visited nodes increases 
significantly. But still, the search performance of BBL is much 
better than BBZ as the lower bound provided by the solution 
of the convex optimization is closer to the true values of the 
subtrees. 

Since calculating lower and upper bounds by means of 
convex relaxation is computationally more demanding than 
calculating the simple bound used for BBZ, the runtime of 
BBZ is lower for short time horizons even if BBZ leads to 
larger decisions trees. But with an increasing length of the 
time horizon, the difference in runtime between BBZ and 
the other BB methods becomes smaller and at some point, 
both methods outperform BBZ. For example, with the current, 
barely optimized implementation based on MATLAB version 
7.9, BBC outperforms BBZ from horizon length A^ = 9 on for 
cost function Ci. It is expected that employing an optimized 
implementation, e.g., with an optimized convex programming 
toolbox like CVX ifTSl . outperforming BBZ occurs for signif- 
icantly shorter time horizons. 

B. Comparison of Objective Function Values 

In Fig. [T] (b), the objective function values of BBC are 
compared with GREEDY, GREEDY*, and CONVEX for the 
costs Ci{N). Both greedy methods are the computationally 
cheapest, but in turn provide highly suboptimal results. Due 
to the myopic planning, GREEDY is not able to anticipate 
the long-term effect of early selecting costly sensors. In this 
simulation example, GREEDY omits measurements at the last 
time steps of the horizon and not in between in order to 
meet the maximum cost constraint. This effect is attenuated 
by GREEDY* thanks to the modified scalar functions gk{ ■ ), 
where the sensor costs are minimized together with the state 
covariance. In doing so, GREEDY* can for example utilize 
sensor 1 instead of the more costly but also more accurate 
sensor 3. This in turn enables GREEDY* to omit less mea- 
surements at the end of the time horizon. 
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CONVEX 
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Fig. 2. Rmse of the target position [x, y]^ wlien GREEDY (red, daslied), 
GREEDY* (blue, dasli-dotted), CONVEX (green, dotted), and BBC (black, 
solid) are employed over a time horizon of length A^ = 10. 

The proposed suboptimal CONVEX method provides sensor 
schedules close to the optimal ones, whereas the computational 
demand is significantly smaller compared to BBC, especially 
for very long time horizons. CONVEX trades scheduling 
quality off against scheduling complexity, which is desirable 
for computationally constrained sensor systems. 

C. Tracking Error 

For a time horizon of length TV = 10 and for the cost 
function Ci{N), 100 Monte Carlo simulation runs are per- 
formed for evaluating the target tracking error when employing 
GREEDY, GREEDY*, CONVEX, and BBC. In Fig. |2] the 
root mean square error (rmse) with respect to the target 
position [x, j/]^ is depicted. Compared to the greedy methods, 
both CONVEX and BBC provide the lowest tracking errors 
together with the lowest uncertainty (consider Fig. [T] (b)), 
whereas BBC performs best. The little bump in the curves 
of CONVEX and BBC around time step fc = 6 results from 
omitting a measurement. Here, GREEDY* has the smallest 
tracking error, which comes at the expense of a higher error 
at the end of the time horizon due to omitting measurements. 

However, in the simulation example considered here, the 
tracking error of GREEDY* is relatively close to the errors 
of the proposed convex optimization based approaches. More 
significant benefits of multi-step sensor scheduling are ex- 
pected for example in scenarios where sensors are temporarily 
unavailable or in scenarios with nonlinear dynamics and sensor 
models (consider for example the results in ifTOl . lfT4l ). The 
application of the proposed approaches to nonlinear models is 
subject of future work. 

VII. Conclusions and Future Work 

Convex optimization is a promising direction for deter- 
mining multi-step sensor schedules. In this paper, a general 
sensor scheduling problem for linear Gaussian systems was 
formulated and the convexity of its relaxation was proven. 
Based on this result, a suboptimal and an optimal scheduling 
approach utilizing convex optimization have been proposed. 



While the suboptimal approach trades estimation quality off 
against computational demand, the optimal one outperforms 
existing optimal algorithms with regard to search speed. 

Compared to existing approaches on sensor scheduling 
via convex optimization, arbitrary linear Gaussian sensor 
scheduling problems are covered. Furthermore, both proposed 
scheduling methods are appropriate for long time horizons and 
many sensors, where choosing the better suited approach for 
a given scheduling problem depends on the requirements on 
estimation quality and computational capabilities. 

Future work is devoted to efficient sensor scheduling for 
very long or even infinite time horizons, where the BB- 
based approach is computationally intractable. Here, model- 
predictive control (also referred to as moving horizon control, 
see for example 1151 ) can be employed. Furthermore, it is 
intended to extend the proposed convex sensor scheduling 
approaches to nonlinear dynamics and sensor models. This 
could be achieved for instance by employing model-predictive 
control in combination with a conversion of the nonlinear 
models into linear ones via linearization, e.g., first-order Taylor 
series expansion or statistical linearization as proposed in llT4l . 
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